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Abstract 

In 1992 Bender, Boettcher, and Lipatov proposed in two papers a new and unusual nonper- 
turbative calculational tool in quantum field theory. The objective was to expand the Green's 
functions of the quantum field theory as Taylor series in powers of the space-time dimension D. 
In particular, the vacuum energy for a massless (j) 2N (N = 1, 2, 3, . . .) quantum field theory was 
studied. The first two Taylor coefficients in this dimensional expansion were calculated exactly 
and a set of graphical rules were devised that could be used to calculate approximately the higher 
coefficients in the series. This approach is mathematically valid and gives accurate results, but 
it has not been actively pursued and investigated. Subsequently, in 1998 Bender and Boettcher 
discovered that "PT-symmetric quantum- mechanical Hamiltonians of the form H = p 2 + x 2 (ix) e , 
where e > 0, have real spectra. These new kinds of complex non-Dirac-Hermitian Hamiltonians 
define physically acceptable quantum-mechanical theories. This result in quantum mechanics sug- 
gests that the corresponding non-Dirac-Hermitian D-dimensional (f) 2 (i<p) € quantum field theories 
might also have real spectra. To examine this hypothesis, we return to the technique devised in 
1992 and in this paper we calculate the first two coefficients in the dimensional expansion of the 
ground-state energy of this complex non-Dirac-Hermitian quantum field theory. We show that to 
first order in this dimensional approximation the ground-state energy is indeed real for e > 0. 




PACS numbers: 11.30.Er, 03.65.-w, 02.30.Fn, ll.10.-z, ll.10.Jj 



I. INTRODUCTION 



Many nonperturbative techniques have been invented to study the properties of quantum 
field theory. In this paper we reconsider a little-known nonperturbative technique that 
was invented in 1992 by Bender, Boettcher, and Lipatov [l], Q] and we use this technique 
to examine the ground-state energy of a massless scalar quantum field theory having a 
self-interaction of the form </> 2 (z</>) € . While the Hamiltonian for this quantum field theory 
is complex and obviously not Dirac Hermitian, the Hamiltonian is 'PT-symmetric. The 
nonperturbative technique developed in Refs. [l[ and (2[ provides support for the conjecture 
that the ground-state energy for this strange looking quantum field theory is real. 

The nonperturbative technique used in this paper is called dimensional expansion and 
the general idea is to expand the Green's functions of a quantum field theory as Taylor 
series in powers of the space-time dimension D. The leading term in such an expansion is 
easy to determine because it requires only that one solve the corresponding zero- dimensional 
version of the quantum field theory; clearly, this is a trivial exercise. The surprising discovery 
that was reported in Refs. [l| and [2[ is that the first- order term (the coefficient of D) in 
the dimensional expansion can also be calculated exactly and in closed form. Furthermore, 
graphical techniques were devised that can be used to obtain approximately the higher-order 
terms in the dimensional expansion. 

In the conventional formulation of quantum mechanics the Hamiltonian H is assumed to 
be Dirac Hermitian, where by this term we mean that H = W and the symbol f represents 
combined matrix transposition and complex conjugation. The condition of Dirac Hermiticity 
allows the matrix elements of H to be complex, but it guarantees that the eigenvalues of H 
are real. However, in 1998 it was shown that it is possible to extend a Hermitian Hamiltonian 
into the complex domain while keeping the eigenvalues of the Hamiltonian real, even though 
the Hamiltonian is no longer Dirac Hermitian 0, H| . These papers showed how to extend 
the harmonic oscillator into the complex domain and thereby obtain the following infinite 
family of complex non- Hermitian Hamiltonians having real eigenvalues: 

H = p 2 + x 2 (ix) e (e>0). (1) 

In Refs. j3[ and 0] it was argued that the reality of the spectra of H in (JT]) was due to the 
unbroken PT symmetry of these Hamiltonians for e > 0. By unbroken PT symmetry, we 
mean that the eigenfunctions of the Hamiltonian are also eigenf unctions of the PT operator. 
The eigenvalues of this class of Hamiltonians are plotted in Fig. [TJ 

Thus, in quantum mechanics the class of physically allowed Hamiltonians can be en- 
larged to include non-Dirac-Hermitian but PT-symmetric Hamiltonians in addition to Dirac- 
Hermitian Hamiltonians, and eigenfunctions (solutions to the Schrodinger equation) are 
treated as functions of complex coordinates [5l-ll0|. Experimental observations of physical 



pnye 

systems described by complex PT-symmetric Hamiltonians are now being reported |lll-ll8 

The fact that PT-symmetric Hamiltonians can define a conventional theory of quantum 
mechanics suggests that it is also possible to extend quantum field theories into the complex 



domain [19[ . The quantum-field-theoretic analogs of the quantum theories described by the 
Hamiltonians in ([1]) are the massless complex scalar quantum field theories whose Euclidean 
Lagrangians have the form 

L = Ud<Pf + (e>0). (2) 



It has not been proved rigorously that the energies of such quantum field theories are real 
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FIG. 1: Energy levels of the Hamiltonians H = p 2 + x 2 (ix) e in ([T]) plotted as functions of the 
real parameter e. Note that there are three regions: When e > 0, the spectrum is discrete, real, 
and positive and the energy levels rise with increasing e. The lower bound of this region, e = 0, 
corresponds to the harmonic oscillator, whose energy levels are E n = 2n + 1. When — 1 < e < 0, 
there are a finite number of real positive eigenvalues and an infinite number of complex conjugate 
pairs of eigenvalues. As e decreases from to —1, the number of real eigenvalues decreases; when 
e < —0.57793, the only real eigenvalue is the ground-state energy. As e approaches — 1 + , the 
ground-state energy diverges. For e < — 1 there are no real eigenvalues. When e > 0, the VT 
symmetry is unbroken, but when e < the VT symmetry is broken. 



but approximate solutions to the Dyson-Schwinger equations indicate that the poles of the 
Green's functions are real 20N24 . 

This paper is organized very simply. In Sec. |H| we review the dimensional-expansion 
techniques that were introduced in Refs. [l[ and [2J and show how to apply them to the 
field theories whose Euclidean Lagrangians are given in fl2]). Then, in Sec. IHII we make some 
concluding remarks and suggest future directions for research. 



II. DIMENSIONAL EXPANSIONS 

Following Refs. [l[ and 0, we derive a formula for the dimensional expansion of the 
logarithmic derivative A(e, D) of the free-energy density F for the Lagrangian L in fl2]): 

dF 

A(e,D) = (2 + e)g—. (3) 

Note that A(e, D), as defined in ([3]), is a dimensionless number for all D and is independent 
of the coupling constant g; A(e, D) is a function only of the dimensionless numbers e and D. 



We seek a dimensional expansion of A(e, D) in the form of a Taylor series in powers of D: 

oo 

A(e,D) = Y,an(e)D n . (4) 



n=0 



A. Zeroth-order calculation 

Calculating the first term in the dimensional expansion (j3J) is easy because it is obtained 
from the zero- dimensional version of the field theory: 

a (e) = A(e,0) 

= —(2 + e)<?-j- log < N / dx exp [— gx 2 (zx) e ] 
^fi 1 I Jc 



(2 + e)<? 



f c dx gx 2 {ix) e exp [— gx 2 (ix) 



f c dx exp [— gx 2 (io 



(5) 



Thus, the first term in the dimensional expansion (j3J) is simply unity. 

It is important to point out that the contour C for the integrals in (J5J) terminates in two 
Stokes' wedges in the complex-x plane. As e increases from 0, these Stokes' wedges (inside 
of which the integral converges) rotate downward and become thinner. At e = the Stokes' 
wedges are centered about the positive-real and negative-real axes and they have angular 
opening 90°. When e > 0, the Stokes' wedges are centered about the two angles 

7re 7re 

bright = -J^Y e aIld #left = -7T + J^^- (6) 

Thus, the two Stokes' wedges are symmetric about the imaginary-x axis. The opening angles 
of the Stokes wedges are vr/(2 + e). Thus, at e = 1 (an i(f> 3 theory) the Stokes' wedges are 
adjacent to and below the real axes and have angular opening 60°. When e > 1, these 
wedges lie below and do not contain the real axis. 

This smooth movement of the Stokes' wedges indicates that the field theories defined by 
the Lagrangians in ([2]) can be viewed as being smooth extensions in the parameter e of the 
Dirac-Hermitian (harmonic-oscillator) Hamiltonian, which lies at e = 0, into the complex 
non-Dirac-Hermitian domain for which e > 0. 



B. First-order calculation 

Next, we show how to calculate the coefficient of D in the dimensional expansion fll]). To 
begin, we add and subtract a mass term in the Euclidean Lagrangian (J2]): 

L = ±(<90) 2 + |m 2 2 + g4> 2 ^Y - im 2 2 , (7) 



where m is a mass parameter. 



The weak-coupling expansion for the Green's functions of this Lagrangian is in principle 
determined by a formal set of Feynman rules: 
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line : 

2 vertex : 
e) vertex : 
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(8) 



Here, the coupling constant g has dimensions of (mass) 2+e De / 2 . Thus, we choose the mass 
parameter in (JTj) equal to the appropriate power of g 



m 
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l/(2+e-£»e/2) 



(9) 



It is important to note that if e is an integer, then there is a diagrammatic expansion defined 
by the rules in (jSJ). However, it is not necessary that e be an integer; it can be any nonnegative 
real number. Thus, unlike the case in Refs. jl[ and 0], the amplitude — i e (2 + e)! g does not 
technically represent a vertex but rather the continuation of a vertex amplitude off integer 
values of e + 2. 

The free-energy density F for the Lagrangian L in flZJ) is given by the standard functional- 
integral representation 



-FV 



exp 



c 



d D xL 



(10) 



where C is the complex contour discussed above, the measure must be properly normalized, 
and V represents the volume of space-time in D- dimensions. We can express the right side 
of ffTUjl as the product of two factors: 



-FV 



J V(f) exp ^~ 


j d D xL^j 
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f c V(j) exp (-/ d D xL) 



(11) 



j c V4> exp (-/ d D xL ) 
Note that Lq is the free (noninteracting) part of L in (JTj): 

L o = i(<90) 2 + imV- (12) 
Next we take the logarithm of (flTT) and express the free energy F as the sum of two terms: 

F = F 1 + F 2 , (13) 

where 



Fx = --log 



V<fr exp 



Uc 



d D xL ( 



and 



Fo = -- 



f c V<f> exp (-/ d D xL) 
j c V<j) exp (-/ d D xL ) 



(14) 



(15) 



The functional-integral expression for F\ is the free-energy density associated with the nonin- 
teracting Lagrangian in (TI2]) . The expression F2 represents the negative sum of all connected 
vacuum graphs except for polygonal loop graphs contributing to the free-energy density of 
the field theory described by the Lagrangian L . This is because the integral in the denom- 
inator of the argument of the logarithm in (I15p eliminates from F 2 precisely those one-loop 



diagrams that contribute to iq. We can in principle calculate these graphs using the Feyn- 
man rules in OH]). We now show how to calculate each of these two quantities. 



Calculation of iq: 

For this calculation we must use the 2 vertex in the Feynman rules in (jSJ). It is possible 
to calculate the free-energy density iq exactly for all dimensions D. 

The free-energy density iq is the negative of the total contribution of all connected 
vacuum graphs. Since these graphs are constructed from a two-point vertex, the graphs are 
all polygons. An n- vertex polygon has symmetry number l/(2n), vertex factor (— m 2 ) n , and 
momentum integral 

w D j tie) 

Hence, apart from an additive divergent zero-point energy contribution that is independent 
of m, we have 

fi = -EsNW/y. (17) 

n=l J y 

Each of the integrals in the summation (1171) is infrared divergent, but if we interchange 
the orders of summation and integration, we see that this divergence evaporates and we 
obtain a finite integral representation for iq: 

F 1 = i(27r)- D /d^log(l + ^). (18) 

This integral converges for < D < 2, but the integral is ultraviolet divergent at D = 2. 

We can evaluate the integral in f fT8|) straightforwardly by introducing polar coordinates. 
The result is 




Thus, taking a logarithmic derivative, we find that the contribution of iq to Ai(e,D) is 
exactly given by 

(27r)^ D / 2 r (1 - \D) . (20) 
The Taylor expansion in powers of D of this quantity begins as follows: 

1_ i[l 0g ( 2vr )- 7 ] J D + 0( J D 2 ), (21) 

where 7 is Euler's constant. These are the first two terms in the dimensional expansion for 
the logarithmic derivative of the free energy of the quantum field theory in ffl2l) . The radius 
of convergence of this dimensional expansion is 2. 

Calculation of F 2 : 

It was shown in Refs. [l[ and that while one might imagine that calculating F 2 is 
difficult because the computation requires the evaluation of the sum of all vacuum graphs, 
the calculation to order D is very easy. The reason for this is that we only need to calculate 
the sum of the vacuum graphs for the case of zero dimensions if we want to know A(e, D) 
to order D. This is true because the contribution of each vacuum graph is a finite quantity 
having dimensions of mP (the dimensions of an energy density) and each graph is finite in 



the limit as D — > 0, unlike the expression for Fi in ( TT9~j) . Thus, each vacuum graph has the 
general form 

m D [constant + 0(D)]. (22) 

Consequently, when we calculate the logarithmic derivative of each graph in F 2 to obtain 
A(e,D), the factor of m D contributes a factor of D. The vacuum graphs need only be 
calculated for the special case D — 0. 

To verify the form in (1221) . we argue as follows: For every Feynman graph we must 
evaluate a D-dimensional momentum integral as required by the Feynman rules in (jSJ). The 
integrand f(p) of this integral is a combination of propagators of the form l/(p 2 + m 2 ). In 
the limit as D — > 0, we may replace f(p) by /(0) because there is no infrared divergence in 
any graph. In fact, the value of the integral as D — >• is just /(0). 

It follows that we can express F2 for small D as a ratio of dimensionless integrals multiplied 
by the mass factor mP: 

^ 2 = ~m D log % Pl ( \ , J - (23) 

We evaluate the integral in the numerator of ( 1231) by reducing it to a real integral along the 
centers of the Stokes' wedges. The result is that 



m D log 



(2/7r) 1/2 r f 1 + 
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2 + ey \4 + 2e 
This gives the order 0(D) contribution from F 2 to A(e,D). 



(24) 



Combination of F\ and F 2 : 

Finally, we combine the contributions to A(e,D) from Fi in (12 ip and from F 2 in (1241) . 
The result is that the dimensional expansion has the form 

A(e, D) = \-\D [4e-^r 2 (l + cos 2 (^)] D + O (D 2 ) , (25) 
and we emphasize that the coefficients of D° and D 1 in this expansion are exact. 



III. FINAL REMARKS AND FUTURE RESEARCH 

In this paper we have used the dimensional-expansion techniques proposed by Bender, 
Boettcher, and Lipatov in 1992 to find exactly the first two terms of the dimensional expan- 
sion (125]) for the ground-state energy density of the Euclidean quantum field theory defined 
by the complex, non-Dirac-Hermitian Lagrangian in The reality of the coefficients in 
this expansion suggests that the entire expansion may well have real coefficients. 

We do not know how to calculate the higher coefficients in the dimensional expansion 
exactly. However, graphical techniques were devised in Refs. [l| and j2[ that can be used to 
find accurate approximations to the higher coefficients. This would be a worthy objective for 
a research project, and an important step in showing that non-Dirac-Hermitian Lagrangians 
can define physically acceptable quantum field theories. 
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